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Abstract 

A polarized gamma ray emission spread over a sufficiently wide energy band from a strongly 
magnetized astrophysical object like gamma ray bursts (GRBs) offers an opportunity to test 
the hypothesis of invisible axion. The axionic induced dichroism of gamma rays at different 
energies should cause a misalignment of the polarization plane for higher energy events relative 
to that one for lower energies events resulting in the loss of statistics needed to form a pattern 
of the polarization signal to be recognized in a detector. According to this, any evidence of 
polarized gamma rays coming from an object with extended magnetic field could be interpreted 
as a constraint on the existence of the invisible axion for a certain parameter range. Based 
on reports of polarized MeV emission detected in several GRBs we derive a constraint on the 
axion-photon coupling. This constraint ^077 

< 2.2 • KT 11 GeV" 1 calculated for the axion mass 
m a = 10~ 3 eV is competitive with the sensitivity of CAST and becomes even stronger for lower 
masses. 
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The Peccei-Quinn (PQ) mechanism [1] remains perhaps the most natural solution to 
the CP problem in QCD. A new chiral Upq(1) symmetry being spontaneously broken at 
some large energy scale, f a , and explicitly broken by the color anomaly at QCD scale 
would allow for the dynamical vanishing of the 9 term and thus the restoration of the CP 
symmetry in strong interactions. The pseudo-scalar field, which drives the relaxation of 
the 9 term to zero is called axion. The most important phenomenological property of this 
axion is its two-photon vertex interaction, which allows for axion to photon conversion in 
the presence of an external electromagnetic and magnetic fields [2] through an interaction 
term ^ 

£<27 = ^da-yy-F^F^ a, g a jjE ■ Ba, (1) 

where a is the axion field, F is the electromagnetic field strength tensor, F its dual, E 
and B the electric and magnetic fields respectively. The axion-photon coupling strength is 
quantified by 

a 1 

where a is the fine-structure constant and £ is an order one 1 parameter which depends on 
the details of the electromagnetic and color anomalies of the axial current associated with 
axion field. For review of the properties of the invisible axion as well as various types of 
constraints on the axion mass and couplings see [3, 4]. 

According to [7] the axion-photon mixing (1) gives rise to vacuum birefringence and 
dichroism, which are qualitatively identical to those arising from the QED magnetized 
vacuum. If one considers a beam of laying polarized monochromatic photons with frequency 
uj and wave vector k propagating in a vacuo along in a uniform magnetic field B laing at a 
nonvanishing angle <fi with the wave vector, due to the birefringence the beam polarization 
becomes elliptical at some distance from the source. On the other hand, the dichroism 
produces a rotation of the ellipse's major axis with respect to the initial polarization. In this 
letter we only consider the axionic dichroism induced rotation angle e of the polarization 
plane of an initially linearly polarized monochromatic beam given in [7, 8]: 

AT 9 2 ail B 2 uj 2 . (m 2 a L\ . 

e = N — Lj -^ sin sin 20, (3) 

ml \ ) 

where m a is the mass of the axion, L is the length of the magnetized region, N is the number 
of passes through the cavity with the magnetic field if for instance a laser experiment like 
PVLAS [9, 10] is considered. The validity of the approximation (3) is provided if the 
oscillation wavenumber 

AL=P^) 2 + ^, (4) 

is dominated by the axion mass term. In fact, (4) pertains to the situation in which the 
beam propagates in a magnetized plasma, which gives rise to an effective photon mass 



1 Two quite distinct invisible axion models, namely the KSVZ [5] (hadronic axion) and the DFSZ [6] 
one, lead to quite similar g ail - 
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set by the plasma frequency cu p \ = yAiran e /m e ~ 3.7 • 10~ n yn e /cm -3 eV, where n e is 
the electron density and m e is the electron mass. Since the dichroism rotation, expected 
to arise in the QED vacuum, and is suppressed with respect to birefringence by factor 
(B/B cr ) 2 [11], where B cr = m e /e ~ 4.4 • 10 13 G, any observable rotation in vacuum should 
be attributed to the axion-photon mixing term (1). 

The polarization of the prompt gamma ray emission has been measured in four bright 
GRBs: GRB021206, GRB930131, GRB960924 and GRB041219a. The first measurements 
made in [12] with Ranaty High Energy Solar Spectrometer Imager (RHESSI) satellite [13], 
found a linear polarization, II = (80 ± 20)%, of the gamma rays from GRB021206 across 
the spectral window 0.15-2 MeV. The analysis techniques have been challenged in [14] and 
defended in [15]. Subsequent analyses made in [16] confirmed the results of [12] but at the 
lower level of significance. Later, in [17] the BATSE instrument on board of the Compton 
Gamma Ray Observatory (CGRO) [18] has been used to measure, for two GRBs, the 
angular distribution of gamma rays back-scattered by the rim of the Earth's atmosphere: 
35% < IT < 100% for GRB930131 and 50% < II < 100% for GRB960924. The analysing 
technique of [17] is only sensitive to the energy range 3-100 keV. Finally, the analysis [19] of 
GRB041219a across the spectral window 100-350 keV has been performed using coincidence 
events in the SPI (spectrometer on board of the INTEGRAL satelite [20]) and IBIS (the 
Imager on Board of the INTEGRAL satelite). The polarization fraction of II = 96±H% 
was determined for this GRB. 

The above mentioned measurements are made using multiple events scattered into 
a detector with geometry distinguishing capabilities (two adjacent detectors in case of 
INTEGRAL or rotating detector in case of RHESSI). Because Compton scatter angle 
depends on the polarization of incoming photons, the time integrated polarization of the 
prompt gamma ray emission in a GRB can be of the order of tens of per cent provided 
that the polarization angle does not vary significantly during the whole duration of the 
GRB across the spectral range analyzed. For example the detected polarization signal in 
RHESSI arises from a correlation between the time dependence of scattered photon flux 
and the angular orientation of the satellite, which varies on a time scale similar to the 
burst duration. The data [21] indicate that a major contribution to the flux comes from 
photons significantly distributed over at least the energy range 0.2-1.3 MeV. 

According to the Hillas [22] 2 diagram showing size and magnetic field strengths of 
different astrophysical object the typical magnetic field in a GRB's engine can be estimated 
as B ~ 10 9 G over a region Lgrb — 10 9 cm. Therefore, one can observe that the constraint 
arises from the fact that if the axionic dichroism induced angle of polarization rotation (3) 
in the given magnetic field were to differ by more then it/2 over the energy range 0.2- 
1.3 MeV, as in the case of GRB021206, the instantaneous polarization in the detector 
would fluctuate significantly for the net time averaged polarization of the signal to be 
suppressed. Similar argument can be applied to the polarization measurements based on 
INTEGRAL and BATSE events. 

To evaluate the bound we assume that in average the magnetic field is sufficiently 

2 See Fig. 10. of [23] for the graphical compilation of [22]. 
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misaligned with the direction of the beam of gamma rays over the magnetized region 
implying that sin 20 ~ 1 and the axion mass is large enough that for a given energy of 
gamma rays uj 

J GRB 

~ 2 

Therefore, the length 



™>I l grb > vr ^ 



can be interpreted as an analog of the passage length akin to that one which defines the 
dimension of a magnetized cavity in a laser PVLAS experiment. Thus the number of 
passages in (3) is given by 

N = (7) 

and thereby when the condition (5) holds the accumulated polarization rotation angle can 
be expressed from (3) as 

" <-ujB 2 . (8) 



/-chi; /A"/-.-. , r y> 
2-7T m 2 , 



Hence, the relative misalignment between the polarization planes of gamma radiation at 
two different energies u)\ and u>2 is given by 

Ae = ^ 9 ^AuB\ (9) 
2tt mi 

where Au = \u>2 — <^i|- The statistical pattern of the time integrated polarization signal 
from a GRB in a detector is preserved for the energy range between i>j\ and c^2 provided that 
the relative misalignment angle (9) is less than 7r/2. This condition can be transformed 
into the bound on the axion-photon coupling as 

m a 

g ail < 7r =. (10) 

By/ AuL GRB 

From the last inequality using B = 1.95 • 10 7 eV 2 and AujLqrb = 5 ■ 10 19 for the given 
magnetic field, energy difference Auj = 1 MeV and the length of magnetized region Lqrb = 
10 9 cm one obtains 3 

771 

g ail < 2.2 • lO" 8 ^ (GeV)" 1 , (11) 

where the inner part of the spectral window 0.2-1.3 MeV {Auo m lMeV) reported in 
polarization analysis of GRB021206 has been used. However, for the axion mass 



™. < (12) 

V L GRB 



3 The values are given in units c = % = 1, so that IT = 195 eV 2 and GeV • cm = 5 • 10 
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Figure 1: The plot of the regions of (m a ,g a77 ) space ruled out by various solar axion 
searches taken from [4] with the bound of the present letter, estimated for the inner part 
of the energy range 0.2-1.3 MeV applied for the polarization measurements of GRB021206 
(dashed dotted line), superimposed. 

the passage length (6) exceeds Lqrb implying that the rotation angle should be expressed 
as a constant: 

e= 9 -fTL(BL GRB ) 2 . (13) 
lb 

The upper edge u>2 ~ 1.3 MeV of the energy range considered, together with the condi- 
tion (12), defines the critical mass 

m crl « 3.5 • 1(T 4 eV. (14) 

Below this mass the rotation angle of the higher energy photons does not depend either 
on their energy or the axion mass and is given by (13), while for the lower energy edge 
of the polarized photons uj\ ~ 0.2MeV the equation (8) is still valid. Therefore, for the 
axion mass bellow the critical one, m a < m crl , the polarization planes misalignment angle 
should be calculated as 

A, = 5^L GRB ^_-_j. (15) 
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The expression (15) holds to be positive down to the mass 



m cr2 = d - U T l » 8 • 10- 5 eV. (16) 
V ^L GRB 

Requiring again that the misalignment angle (15) does not exceed 7r/2 in the axion mass 
range between m cr i and m cr2 one arrives to a bound, which can be well approximated by 
a constant 



9ayy < « 5 • lO" 12 (GeV)- 1 . (17) 

For the masses m a > 8 • 1CT 5 eV the rotation of the polarization plane is given by (13) 
and misalignment does not appear making the bound undefined for the axion mass bellow 
m cr2 . 

In Fig. 1. we show the bounds (11) and (17) superimposed on the recent results of 
CAST [24] and other axion helioscope experiments [25, 26, 27, 28, 29]. One can see that 
the bound obtained from the lack of a substantial misalignment of polarization planes of 
gamma radiation at different energies within the energy range reported by polarization 
measurements of GRB021206 is competitive with the current sensitivity of CAST for the 
axion masses below 10 -3 eV. Moreover, one observes that it seems to be very unlikely that 
a laser experiment like PVLAS could find a signal [9] corresponding to g ail ~ 10" 5 GeV -1 
for the meV axions. Indeed, the lack of the observation of the previously claimed rotation 
signal [9] has been established recently by the PVLAS team after a substantial upgrade of 
the facility has been made [10]. Also the PVLAS anomaly has not found a support in the 
results of a pulsed "light shining through a wall" experiment [30]. 

Since according to [31] the electron number density in a GRB's environment can be 
estimated as n e ~ 10 10 cm~ 3 the expression (4) is still axion mass dominated down to 
m a pa m CT for the given magnetic field, the energy of the gamma radiation, uo ~ 1 MeV, 
and constraints on g ail calculated from (11) and (17) . Therefore, the validity of the 
approximation (3) holds in the range of the parameters the bound in Fig.l. is imposed. 
The minimal time scale of variability of GRBs light curves is estimated to be about 0.1 
sec 4 . This implies that the typical extention of the GRB's engine is indeed compatible 
with Lqrb ~ 10 9 cm we used for the evaluation of the bound. Conservation of magnetic 
field energy at the rest wind frame of fireball shell model of the GRB's engine [31] implies 
at any radial distance r, in the fireball environment, Airr^B 2 = 4nr 2 B 2 , leading to the 
relation B = B (r /r), where B and r are the magnetic field strength and the size of the 
central part of the fireball. Typically the central part of the fireball can be represented by 
a neutron star of radius r ~ 10 6 cm with magnetic field of B pa 10 12 G. Therefore the 
strength of the magnetic field at the distance r = Lqrb corresponds to B pa 10 9 G, which 
is in a good agreement with the values taken from [22, 23]. 



The limit obtained becomes by factor ylMeV/Au;/^ weaker if we apply the width 
Acjj pa 250 keV of the energy bands for GRB041219a detected by INTEGRAL or Au B pa 



'See, for example, the analysis in [32]. 
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100 keV for GRB930131 and GRB960924 detected by BATSE. This implies that g ail < 
4.4- 10~ n GeV- 1 and g ail < 6.9- 10~ n GeV" 1 for INTEGRAL and BATSE measurements 
respectively calculated for the axion mass m a = 10~ 3 eV. A stronger constraint could be 
obtained from the GRBs considered by taking into account more precisely the spectral 
characteristics of the signal and statistical criteria for loosing the polarization pattern in 
the detectors. Of course, the robustness should find its confirmation in further detection of 
gamma polarized signals from other GRBs in the similar energy range. In principle, apart 
from the mentioned instruments SWIFT [33] satellite is also capable of polarimetry [34] in 
300 keV-10 MeV energy band. 

In recent years, numerous efforts have been initiated to develop instruments with the 
sensitivity required for astrophysical polarimetry over 100 eV to 10 GeV band [35]. Time 
projection chambers (TPCs), with their high-resolution event imaging capability, are an 
integral part of some of this efforts. At the energy band 300 keV-10 MeV, Compton 
polarimeters based on the use of high-Z scattering elements (coupled with high-Z absorbers) 
become viable. For example, the Ge double scatter approach used by RHESSI becomes 
most effective at energies above 300 keV. Liquid rare-gas TPCs are being pursued as 
large effective area Compton telescopes, while electron-tracking gas TPCs are a component 
of Compton telescopes with the lowest background, especially for the most polarization 
sensitive events. In the pair production regime, which is effective in the energy range 
2 MeV-10 GeV, only TPCs currently offer the hope [35] of tracking electron-positron pairs 
or recoil electrons with the accuracy and efficiency required for astronomical polarimetry. 
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